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We have systematically studied superconducting instabilities in the repulsive Hubbard model 
for d-wave and p-wave pairing in various 2D and 3D lattices. Using fluctuation exchange ap- 
proximation, we consider 3D face centered cubic lattice, 3D body centered cubic lattice, 3D 
simple cubic lattice, 2D square lattice and 2D triangular lattice, where either strong ferro- or 
antiferromagnetic spin fluctuation is present. We show that (i) d-wave instability mediated by 
antiferromagnetic spin fluctuations is stronger than p-wave instability mediated by ferromagnetic 
spin fluctuations both in 2D and 3D, and (ii)d-wave instability in 2D is much stronger than that 
in 3D. These amount that the "best" situation is the antiferromagnetic- fluctuation mediated in 
2D as far as the single-band Hubbard model on ordinary lattices are concerned. 

KEYWORDS: Hubbard model, anisotropic superconductivity, fluctuation exchange approximation 
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§1. Introduction 

There has been an increasing fascination with elec- 
tronic mechanisms for pairing since the high=T c super- 
conductivity was discovered in the cuprates.cP Among 
various mechanisms, the pairi ng mediated by spin fluc- 
tuations has been proposeafiu DLmLI and have been ex- 
amined intensively. 

The pairing interaction, denoted by V^, resulting 
from exchanging spin fluctuations is dominated by the 
spin susceptibility which has a strong wavenumber de- 
pendence, so that in turn has a strong wavenumber 
dependence. This determines the symmetry of the pair- 
ing, since this dictates the dominant pair-hopping pro- 
cesses on the Fermi surface. There, the dominant mode 
of spin fluctuations, hence the symmetry of the pairing, 
must be sensitively correlated with the lattice structure 
and filling in general. For example, the respulsive elec- 
tron correlation causes the face centered cubic (FCC) 
lattice with low density of electrons to have strong ferro- 
magnetic fluctuations, or the half- filled body centerd cu- 
bic (BCC) lattice to have strong antiferromagnetic fluc- 
tuations as recently confirmed]!!) Thus, it is an interest- 
ing question to ask how the superconductivity would be 
correlated with the lattice structures. 

As for the singlet, anisotropic {e.g., d-wave) pairing in 
2D systems, calculations on a microscopic level have been 
performed: For the repulsive Hubbard model, which is a 
simplest possible model of the electron correlation, the 
fluctuation exchange, approximation (FLEX) developed 
by Bickers et al.Ui3 has been employed to show the 
occurrence. of the d-wave superconductivity in a square 
latticeOEP As for other 2D lattices, several authors 
have suggested the occurrence of d-wave superconduc- 
tivity in the Hubbard model on 2D anisotropic triangu- 
lar lattice, which-ue^peats-an-organic superconductor 
(BEDT-TTF) 2 XBByil§By> These results indicate 
that T c for the superconductivity near the antiferromag- 



netic instability in 2D is ~ O(0. Oli) (t: transfer integral), 
i.e., two orders of magnitude smaller than the original 
electronic energy, but still 'high T c ' ~ 0(100 K) if we 
take t ~ 0(1 eV) for the cuprates (while T c ~ O(10 K) 
if we take t ~ O(0.1 eV) for the organic^onductors) . 

As for 3D systems, Scalapino et alt3 showed that 
paramagnon exchange near a spin-density wave insta- 
bility gives rise to a strong d-wave pairing interaction for 
the 3D Hubbard model on the simple cubic_Lattice, but 
T c was not discussed there. Nakamura et aJny extended 
Moriya's spin fluctuation theory of superconductivity^ 
to 3D systems, and concluded that T c is similar between 
the 2D and 3D cases provided that common parameter 
values (scaled by the band width) are taken. However, 
the parameters there are phenomelogical ones, so it is 
not clear whether the result remains valid for microscopic 
models. 

The possibility of the triplet pairing, on the other 
hand, has been investigated since the 1960's £pr su- 
perfluid 3 He,clP a heavy fermion ,ss"s±£m UPtapix or 



most recently, an oxide Sr2Ru04Ej , E-3 ) The experi- 
mental results suggesting p-wave pairing in these ma- 
terials have stimulated theoretical studies for electron- 
repulsion originated superconductivity. For the electron 
gas model, Chubukov extended the Kohn-Luttinger the- 
orem,^ which asserts that the repulsively interacting 
electron gas should be instable against pairing forma- 
tions at low enough temperatures, to p-wave pairing 
for 2D and 3D electron gas in the dilute limit bv_an- 
alyzing_the singularity of the scattering amplitude OlT 
TakadaEf discussed the possibility of p-wave supercon- 
ductivity in the dilute electron gas with the Kukkonen- 
Overhauser model.Ef This model considers the effective 
electron-electron interaction composed of the bare inter- 
action and the interactions mediated by charge/spin fluc- 
tuations that contains the so-called local-field correction. 

As for lattice systems, which is the subject of the 
present paper, 2D Hubbard model with large enough 
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next-nearest- neighbor hopping (t') has been shown .-tp 
exhibit p-pairing for dilute band fillings. HlubinatlP 
reached a similar conclusion by evaluating the supercon- 
ducting vertex in a perturbative way.EP However, the 
energy scale of the p-pairing in the Hubbard model, i.e., 
T c , has not been questioned so far. Using a phenomcno- 
logical approach, Monthoux and LonzaricbEJ' have re- 
cently concluded for 2D systems that the d-wave pairing 
is much stronger than p-wave pairing. 

Given these backgrounds, from microscopic view of 
point, we investigate in this paper various 2D and 3D 
lattice structures (square, triangular, simple cubic, BCC, 
FCC) with systematically varied next-nearest neighbor 
hopping to tune the dispersion and varied band filling. 
Specifically, we address the following fundamental ques- 
tions: 

(i) Can the pairing instability in 3D be stronger than 
that in 2D? 

(ii) Can the pairing instability with other symmetry, i.e., 
spin-triplet p-pairing in the presence of ferromagnetic 
spin fluctuations, become strong? 

A p&ct of the present study has been briefly re- 
ported,!^ while here we extend the calculation to var- 
ious cases in order to extensively confirm our previous 
conclusions. Namely, we study possibility of d-wave su- 
perconductivity for nearly half-filled square, simple cubic 
(SC), and BCC lattices where strong antiferromagnetic 
fluctuations are present, along with possibility of p- wave 
superconductivity for low density square lattice with sig- 
nificant next nearest-neighbor hopping, quarter-filled tri- 
angular lattice, and low-density FCC lattice where fer- 
romagnetic spin fluctuations should be dominant. 

We employ the FLEX approximation, which enables 
us to handle strong spin fluctuations. As seen from the 
result, to determine the best situation for superconduc- 
tivity in the Hubbard model we have to consider various 
factors such as the form of the pairing interaction and 
the energy/momentum-dependence of Green's functions 
(lifetime of quasi-particles, etc), so that the way in which 
the pairing instability is correlated with the lattice struc- 
ture is a highly nontrivial problem which is by no means 
predictable from the outset. 

Here we shall show that 

(i) d- wave instability mediated by antiferromagnetic spin 
fluctuations is stronger than p-wave instability mediated 
by ferromagnetic spin fluctuations both in 2D and 3D, 
and 

(ii) pairing instability in 2D is much stronger than that 
in 3D. 

Thus the 'best' situation for the spin fluctuation medi- 
ated pairing is suggested to be the 2D case with dominant 
antiferromagnetic fluctuations as far as the single-band 
Hubbard model on ordinary lattices are concerned. 

§2. Formulation 

2.1 Model Hamiltonian 

We consider the single-band Hubbard model, 

U = *«' c L c .jV + U Y1 ™it n ii> (!) 



where c ia creates an electron at the i-th site with spin 
cr, ilia = c\ a Ci a is the number operator. We consider 
the transfer between second-nearest neighbors, i,j = t' , 
along with t(= 1 hereafter) for nearest neighbors. 

The energy dispersions for square and triangular lat- 
tices are, 

2 

e D (k) = 2i^cos(fc,) + 4t' cos(/ci) cos(/c 2 ), (2) 
i=l 
2 

e A (k) = 2i^cos(fc l ) + 2icos(£;i + k 2 ), (3) 

i=l 

respectively, where we take a square Brillouin zone for 
the latter case as well by inserting diagonal transfers in 
the square lattice. 

The dispersions for SC, BCC and FCC lattices are 
given as 

3 

e sc (k) = 2< cos hi + 4t' ^ cos fe * cos k 3 > ( 4 ) 

i— 1 i<j 

3 

e BCC (k) = 8£cos£;i cosfecosfcs + 2t' Y"cos(2fci), (5) 

i=l 

3 

e FCC (k) =4i^cosfciCosfcj + 2i'^cos(2/s i ), (6) 

i<j i=l 

respectively, where (ki,k2,ks) = {k x ,k y ,k z ). To facili- 
tate the calculation here we take a cubic Brillouin zone 
(— 7r < ki < 7r) by considering a four(two) equivalent, 
interpenetrating BCC(FCC) lattices. 

The density of states for the non-interacting case is dis- 
played in Fig. 1 for the square lattice with t' = 0.0,0.5, 
triangular lattice, SC lattice with t' = 0.0,0.1,0.3, FCC 
lattice with t' = 0.0,0.5 and BCC lattice with t' = 
0.0,0.1. 



Fig. 1. The density of non-interacting states for (a) square lattice 
with t' = 0.0 (solid line), t' = 0.5 (dashed line), (b) triangular 
lattice, (c) SC lattice t' = 0.0 (solid line), t' = 0.1 (dashed line), 
t' = 0.3 (long dashed line), (d) BCC lattice with t' = 0.0 (solid 
line), t' = 0.1 (dashed line), and (e) FCC lattice with t' = 0.0 
(solid line), t' = 0.5 (dashed line). 



2.2 Method 

The FLEX, introduced by Bickers et al.pt$ starts 
from a set of skeleton diagrams for the thermodynamic 
potential, <&, introduced by Luttinger and Ward. $ is a 
functional of Green's function, 67, and a (/c-dependent) 
self energy can be computed by a functional derivative 
of $ with 67. r Ehr_.the calculation of <f> the idea of Baym 
and KadanoffESLj) of taking an important series of dia- 
grams is employed. Hence the FLEX approximation is a 
self-consistent perturbation approximation with respect 
to on-site interaction U. The FLEX approximation is 
suitable for the analysis of Fermi liquid with strong spin 
fluctuations. 
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The self energy is obtained as 

s ( fc ) = ^E G ( fc -<^ (1) (<z), ( ? ) 

in which RPA-type bubble and ladder diagrams are col- 
lected for the interaction, 



for the singlet pairing and 



V {1 \q)= l -U 2 XiA<l) 



1 



+ 2 C/2 Xirr(<7) 



1 + UXirM) 
1 



1 - ^Xirr(g) 



-U 2 X ilr (q) 



with 



Xirr(g) 



Here we have denoted q = (q, ie u ) and k = (k, iuj n ), 
e v = 2ttvT is the Matsubara frequency for bosons while 
uj n = (2n— 1)ttT for fermions, and N is the total number 
of sites. For simplicity, we neglect the diagrams in the 
particle-particle channel. 

The Dyson equation is written as 

G(k, u n y l = G°(k, ujnf 1 - £(k, Un ), (8) 

where G° is the bare Green's function, 

G°(^n) = T- / , (9) 

\uj n + fj, — e k 

with e£ being the energy of a free electron. We have 
solved the equations (7) ~ (||) by setting the chemical 
potential /i so as to fix the density of electrons. 

To obtain T Cl we solve, with the power method, the 
eigenvalue (Eliashberg) equation, 

AS(2) ( fc ) = §E s(2) ( fc ')i G ( fc, )i 2 ^ (2) (fc,fc'), (io) 

k' 

where G(k) the dressed Green's function, and S^ 2 - 1 (k) the 
anomalous self energy, and T = T c corresponds to the 
point at which the maximum eigenvalue A]yiax reaches 
unity. 

The interaction originates from the transverse 

spin fluctuations, longitudinal spin fluctuations and 
charge fluctuations, namely, 



V^(k,k') 



-U 2 
1 



TjXchO - k') 



X zz (k-k')+ X ± (k + k') 



U^xlAk - k') 



l-U 2 xUk~k>) 



U 2 X ur(k + k') 



V (2 Hq) = \ 



u 2 xM 



i - UxUq) 



u 2 x- m (g) 

1 + ^Xirr(g) 



(12) 



1 - U X ur(k + k') 



where Xch is the charge susceptibility, while X zz (x ± ) is 
the longitudinal (transverse) spin susceptibility in the 
RPA form where the dressed Green's function is used. 

Since we have T,^(k) = 'S^(-k) for the spin-singlet 
pairing, whereas T,^(k) — — £( 2 )(— k) for the spin- 
triplet pairing, V^- 2 \k, k') becomes a function of k — k' = 
q with 

U 2 X irr(q) 1 , 1 T U 2 XuM 



1 - UxiAq) 



l + U X irr(q) 



(11) 



for the triplet pairing. We take N — 64 2 sites with n c — 
2048 Matsubara frequencies for 2D square lattice and 
with n c = 1024 Matsubara frequencies for 2D triangular 
lattice, and N = 32 3 with n c = 1024 for 3D. 

§3. Results 

3.1 Square lattice with antiferromagnetic spin fluctua- 
tions 

Let us start with the case of nearly half-filled square 
lattice, which has strong antiferromagnetic fluctuations. 
In this case, T c ~ 0.02 for (i-wave pairing should be ob- 
tained within the FLEX approximation as mentioned in 
the Introduction. We first present the result for t' = 0.0, 
n = 0.85 (0.15 holes doped) and U = 4, which will serve 
as a reference for other lattices. 

In Fig. 2, we show |G(k, i7rfceT)| 2 , a quantity which 
appears in the right-hand side of the Eliashberg equation 
(10). We can see that |G| 2 takes large values (~ 8.0) 
around the Fermi surface. Note that the smaller the self 
energy correction, the peak of |G(k, iirkBT)\ 2 becomes 
larger and the Fermi-liquid picture becomes more valid. 
Also, a large |G(k, iirkBT)\ 2 favors superconductivity as 
we can see from the Eliashberg equation (10). 

In Fig. 3, we plot the susceptibility in the RPA form, 
x(k, 0) = Xirr/(1~ Uxirr), as a function of the wave num- 
ber for T = 0.03. Dominant antiferromagnetic spin fluc- 
tuations are seen as x peaked around k = (ir, tt). A large 
peak in x(k, 0) should imply a large pairing interaction 
as seen in the Eliashberg equation (10). The spread of 
x(k, 0) around the peak in the momentum space is also 
important as we shall come back later, because it mea- 
sures the fraction of effective channels that contribute to 
in the Eliashberg equation (10). 

We also plot in Fig. 3 Imx(kMax,^) against u>, where 
kMax is the momentum for which x(k, 0) becomes max- 
imum and we have normalized Imx(kMax,w) with the 
maximum value. We have obtained the dependence on 
real u with an analytic continuation in the Pade approxi- 
mation. cIP The spread of Imx(kMaxj UJ ) around the peak 
in the ui sector is an electronic couterpart to the Debye 
frequency lob in the BCS theory for the electron-phonon 
system, so this quantity is another important factor in 
the pairing instability. 

In Fig. 4, we plot Aiviax as a function of temperature T 
(normalized by t) along with the reciprocal of the peak 
value of x(k, 0). The pairing instability is measured by 
how Am & x is close to unity, while 1/x — > signifies the 
magnetic ordering. AMax is seen to be extrapolated to 
unity at T ~ 0.02, in accord with previous results. 13 

Here let us comment on a finite T c in 2D systems. As 
usually done, it is taken to be a measure of T c when the 
layers are stacked with the Josephson coupling. How- 
ever, if we take account of the superconducting fluctu- 
ations rigorously, T c for a purely 2I^s^teni must be 
zero according to Mermin's theoremSSEf To judge 
whether a mean field treatment of pairing instabilities 
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is adequate will require an evaluation of the coherence 
length of the pairing. Although this has been done for 
the conventional phonon mechanism of superconductiv- 
ity with Gor'kov's argumentJlJ' such an evaluation is 
rather difficult in the present case of an electron mech- 
anism, since the pairing potential strongly depends on 
frequency and wave-number. 

We can on the other hand study the effect of a weak 
three dimensionality on the AMax within the present for- 
malism. Specifically, we introduce an inter-layer hop- 
ping, t z . Figure 5 shows the results for AMax for a 3D 
anisotropic cubic lattice with the inter-layer hopping var- 
ied over 0.0 < t z < 1.0 on top of t x = t y = 1.0. We can 
see that weak but finite t z < 0.3 does not appreciably 
change the behavior of AMax- Thus, provided the coher- 
ence length in z-direction is large enough for t z ~ 0.3 to 
validate the mean-field treatment, T c obtained here for 
2D can indeed be used as a measure of T c when a weak 
three dimensionality is present. 

Fig. 2. The squared absolute value of the Green function for the 
smallest Matsubara frequency, iw n = ink^T, against wave num- 
ber for the Hubbard model on a square lattice with t' = 0, 
n = 0.85 and (7 = 4. 



Fig. 3. (a)Imx(kjv[axi w ) (normalized by its maximum value) as 
a function of ui/t for the Hubbard model on a square lattice with 
t' = 0.0, n = 0.85, U = 4 and T = 0.03. (b)x(k, 0) as a function 
of wave number for the Hubbard model on a square lattice with 
t' = 0.0, n = 0.85, T = 0.03 and [7 = 4. 



Fig. 4. The maximum eigenvalue of the Eliashberg equation 
(solid lines) and the reciprocal of the value of \ (dashed lines) at 
the antiferromagnetic peak against temperature for the Hubbard 
model on a square lattice with t' = 0, n = 0.85 and U = 4. 



Fig. 5. The maximum eigenvalue of the Eliashberg equation 
against temperature for the Hubbard model on an anisotropic 
cubic lattice with t x = t y = 1 and t z = 0.0 ~ 1.0, n = 0.85 and 
[7 = 4. 



3.2 t-t' square lattice with ferromagnetic spin fluctua- 
tions 

We move on to a 2D case with dominant ferromagnetic 
spin fluctuations, where triplet pairing is expected. The 
situation for which the ferromagnetic fluctuations be- 
come dominant has extensively been investigated for the 
Hubbard model with various approaches, and one guid- 
ing principle is that a large density of states at the Fermi 
level located near the bottom of the band should favor 
ferromagnetism for sufficiently strong electron-electron 



repulsion. 

For the square lattice, such a situation may be realized 
for relatively large 0.5) for dilute electron densities. 
In this case, the divergence (van Hove singularity) in the 
density of states, which is at the center and has a func- 
tional form |lni?| for t' — 0, shifts toward the band bot- 
tom with D(E) ~ 1/^/(2 + E)\\n{2 + E) | for t' = 0.5, so 
that the density of states at the Fermi level beepmes large 
for the dilute case, favoring ferromagnetism.cif It has 
in fact been shown from quantum Monte Carlo studies 
that the ground state is fully spin-polarized for t' = 0.47, 
n ~ 0.4.c3'E!P So we explore the possibility of p-wave in- 
stability associated with this ferromagnetism, and com- 
pare the result with that of d-wave instability associated 
with the antiferromagnetism discussed in the previous 
subsection. 

First, in Fig. 6, we show AMax f° r the density varied 
over 0.3 < n < 0.6 and t' varied over 0.3 < t' < 0.6 for 
U = 4, 6 with T ~ 0.03. We can see that Am & x becomes 
largest for n — 0.3 and t' — 0.5 for both U = 4 and 6, so 
we take these parameter sets. 

In Fig. 7, we plot 0) for U = 4 as a function 
of momentum along with Imx(kMax,^) as a function of 
ui. The peak of x(k, 0) is indeed located at the T point 
(k = (0,0)) indicating ferromagnetic fluctuations. The 
frequency spread is similar to the case of the nearly half- 
filled t' = square lattice. In Fig. 8, we plot Am & x as a 
function of T. We can see that AMax is much smaller 
than that in the antiferromagnetic case, Fig. 4. 

Fig. 6. The maximum eigenvalue of the Eliashberg equation, 
Ajjax, for a square lattice with the density varied over 0.3 < 
n < 0.6 and t' varied over 0.3 < t' < 0.6 for U = 4, 6 with 
T = 0.03. 



Fig. 7. (a)Imx(kjviaxi (normalized by its maximum value) as 
a function of ui/t for the Hubbard model on a square lattice with 
t' = 0.5, n = 0.3, [7 = 4 and T = 0.03. (b)x(k, 0) as a function 
of wave number for the Hubbard model on a square lattice with 
t' = 0.5, n = 0.3, T = 0.03 and U = 4. 



Fig. 8. A plot similar to Fig. 4 for the Hubbard model on a 
square lattice with t' = 0.5, n = 0.3 and [7 = 4. The dashed line 
denotes the reciprocal of the ferromagnetic peak in x- 



Fig. 9. A plot similar to Fig. 2 for the Hubbard model on a 
square lattice with t' = 0.5, n = 0.3 and [7 = 4. 



3.3 Why is p-pairing weaker than d-pairing? 

The present result that the p-pairing has a lower T c 
contrasts with a naive expectation from the BCS picture, 
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in which the T c should be high for a large density of states 
at the Fermi level. We may trace back the reason why 
this does not apply as follows. 

First, if we look at the dominant (oc 1/[1 — U\i II (q)]) 
term in the pairing potential itself in eqs. (11) and 
(12), the triplet pairing interaction is only one-third of 
the singlet pairing interaction. So this should be one 
reason. 

On the other hand, the large density of states, namely 
the flatness of the band around the Fermi level, is re- 
flected to the fact that |G| 2 (Fig. 9) forms an almost flat 
plateau in a large portion of the Brillouin zone. If the 
maximum value of \G\ 2 is large enough to compensate the 
disadvantage of the one-third , a large A may emerge. 
However, we can see that \G\ 2 is much smaller than that 
in the antiferromagnetic case (Fig. 2), which implies that 
the self energy correction is large. Even when we take 
a larger repulsion U to increase the triplet pairing at- 
traction (i.e., to increase the susceptibility), this makes 
the self-energy correction even greater, resulting in only 
a slight change in A. 
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refJ^ that closer the peak of the density of states is to 
the band bottom, the stronger the ferromagnetic fluctu- 
ations tend to be. 

Despite this difference in ferromagnetic fluctuation, 
the AMax in the present case is comparable to that in 
the square lattice with t' . This may be because the fre- 
quency spread in the present case is quite large as seen 
in Fig.lO(a). 

Fig. 10. (a)Imx(kMaxi w ) (normalized by its maximum value) as 
a function of uj/t for the Hubbard model on a triangular lattice 
with n = 0.5, U = 8 and T = 0.03. (b)x(k, 0) as a function of 
wave number for the Hubbard model on a triangular lattice with 
n = 0.5, T = 0.03 and U = 8. 



Fig. 11. A plot similar to Fig. 4 for the Hubbard model on a 
triangular lattice, with n = 0.5 and U = 8. The dashed line 
denotes the reciprocal of the ferromagnetic peak of \. 



3.4 2D triangular lattice 

Next, we discuss the case of 2D triangular lattice. As 
mentioned in the Introduction, the half-filled case with 
dominant antiferromagnetic fluctuaticpSp-h^p^keady 
been discussed by a number of authors. E3lI5II§ll30E§) 

Here, we focus on the quarter-filled isotropic triangu- 
lar lattice, where we can expect ferromagnetic fluctua- 
tions, and hence p-wave pairing. As we can see in Fig. 
1, the density of states for the triangular lattice has a 
sharp peak (with D(E) ~ \ln(E + 2)|) near the bottom 
of the band, so a dominant ferromagnetic fluctuation is 
expected if the Fermi level is located at the peak. Hanish 
et als3 studied the instability of the fully-ferromagnetic 
state of Hubbard model on the triangular lattice for large 
U and concluded that the ferromagnetic state is stable 
for n ~ 0.5 (quarter- filled) . 



In this situation, we have calculated A 



Mm 



for T 



0.03, U — 4, 8, 12. Since Am &x takes similar values for 
U = 4, 8, 12, we take U = 8. 

In Fig. 10, we plot the wave-number dependence of 
x(k) and frequency dependence of Imx(kMax,^)- We 
can see that the peak of x(k, 0) is located around the T 
point. In Fig. 11, we plot AMax along with the recipro- 
cal of the peak value of x(k, 0) as a function of T. We 
can see that the p-wave instability here is again much 
weaker than the d-wave instability for the nearly half- 
filled square lattice (Fig. 4). The one-third factor dis- 
cussed in the previous subsection should again be one 
factor. The second factor is also involved, namely, the 
maximum value of U 2 \G(k, ik^T)^ is small (~ 1.0 x 8 2 , 
see Fig. 12) compared to that in Fig. 2 (~ 8.0 x 4 2 ), so 
the quasi-particles are short-lived. 

We note that the ferromagnetic fluctuation in the tri- 
angular case is much weaker than in the square lattice 
with large t' (compare Figs. 10(b) and 7(b)), although 
we have taken a larger U here. This may be because the 
peak in the density of states in the triangular lattice is 
situated slightly above the bottom (compare Figs.l (a) 
and (b)). Onoda and present authors have suggested in 



Fig. 12. A plot similar to Fig. 2 for the 2D Hubbard model on 
triangular lattice with n = 0.5, and U = 8. 



3.5 3D simple cubic lattice 

Let us now move on to the case of 3D systems, we 
first discuss the possibility of d-wave pairing in the nearly 
half-filled Hubbard model on the 3D SC lattice. In this 



case, we find that the representation 



V etc, 



with the gap function A(k) ~ cosfc x — cosfc y etc) of Oh 
grouped has the largest AMax, so we concentrate on this 
pairing symmetry. 

In Fig. 13, we show Am & x for the density varied over 
0.75 < n < 0.9 and t' varied over -0.5 < t' < +0.4 for 
U = 4,6, 8, 10, 12 with T = 0.03. We can see that among 
these parameter sets, AMax becomes largest for n = 0.8, 

t! = -0.2 0.3 and U = 8 ~ 10, so we take these 

parameter sets. 

In Fig. 14(a), we show the momentum dependence of 
X(k,0) for U = 8, n = 0.8, t' = -0.3, T = 0.03. We 
can see that the peak of x(k, 0) is located around the K 
point (tt, it, 7r) as expected for the antiferromagnetism. 

In Fig. 15, we plot AMax along with the reciprocal 
of the peak value of x(k, 0) as a function of T for t' — 
—0.2, —0.3 ,U — % and n = 0.8. We see that Am 3X does 
not become very close to unity in the range calculated 
here. For T < 0.02, the result obtained for N = 32 3 and 
n c = 1024 is not convergent enough with respect to the 
system size and the number of Matsubara frequencies, 
so we show the result for a larger n c — 2048 in the inset, 
which suggests that Am & x tends to increase with n c . We 
have also performed a calculation for N — 16 3 , and found 
that AMax also increases with N, so a finite T c (< 0.01) 
may be obtained at least for t' = —0.3, U — 8, n = 
0.8 in the limit of large N and n c . At any event, T c 
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is significantly smaller than in the square lattice with 
strong antiferromagnetic fluctuations (Fig. 4). 

Fig. 13. A plot similar to Fig. 6 for an SC lattice, with the 
density varied over 0.75 < n < 0.9 and t' varied over —0.5 < 
t' < 0.4 for U = 4,6,8, 10, 12 with T = 0.03. Where the curves 
are truncated the system becomes antiferromagnetic (i.e., Uxiir 
becomes unity). 



Fig. 14. (a) 0) as a function of wave number for the Hubbard 
model on an SC lattice with t' = -0.2, n = 0.8, T = 0.03 and 
U = 8. (b)Imx(kMaxi w ) (normalized by its maximum value) as 
a function ofui/t for the Hubbard model on an SC lattice(dashcd 
line) with t' = —0.2, n = 0.8, U = 8, as compared with the 
Hubbard model on a square latticc(solid line) with tl = 0, n = 
0.85, U = 4 for T = 0.03 (left frame, where the 2D and 3D 
results almost overlap) or T = 0.04 (right frame). 



Fig. 15. A plot similar to Fig. 4 for the Hubbard model on an SC 
lattice with t' = -0.2, -0.3, n = 0.8 and [7 = 8. The inset shows 
the results for a larger number (=2048) of Matsubara frequencies 
for *' = -0.3. 



Fig. 16. A plot for the squared Green function against k x and k y 
with k z = 0, it/2, it for the Hubbard model on an SC lattice with 
t' = -0.2, n = 0.8 U = 8, T = 0.03. 



3.6 Why is d-pairing stronger in 2D than in 3D? 

So the d-wave instability is decidedly stronger in 2D 
than in 3D as far as the square and simple cubic lattices 
are concernded, and the question is: what are physical 
reasons for that. We can pinpoint the origin by looking 
at the factors involved in the Eliashberg equation (10), 
i.e., (a) the factor [/ 2 |G| 2 , (b) the summation over the 
frequency, and (c) the summation over the momentum. 
In addition, the factor is of course important in 

the equation, but in the following we compare 2D and 
3D in the situation where the maximum value of x(k, 0) 
(that determines V^) is similar between the two cases 
to concentrate on the factors (a)(b)(c). 

In Fig. 16 we plot |G| 2 in 3D for k z = 0,7r/2,7r as 
a function of k x and k y for U = 8, n = 0.8, t' = —0.3, 
T = 0.03. We can see that the maximum value of U 2 \G\ 2 
is greater in 3D than in 2D. 

Fig. 14(b) displays Imx(k Max ,o;) as a function of uj. 
The figure compares the SC lattice (t' = —0.2, n = 
0.8, U = 8) with a typical square lattice with t' = 0, 
n = 0.85 and U — 4 having a similar magnitude of x 
at T = 0.03. We can see that the frequency spread of 
Im%(kMax,w) is similar between 3D and 2D, so the fac- 
tors (a)(b) can be excluded from the reason for the 2D- 
3D difference. Note that if the frequency spread of the 



susceptibility is scaled not by t but by the band width, 
as Nakamura et aE$ have assumed, AMax would have 
become larger. 

If we turn to the momentum sector in the suscepti- 
bility, x(k, 0), Fig. 14(a) shows that the width, a, of 
the x(k, 0) peak in each momentum component is sim- 
ilar to those in 2D displayed in Fig. 3, where the main 
contribution of to A is confined around (n, it) in 
2D or (-7T, 7T, 7r) in 3D. Since the right-hand side of the 
Eliashberg equation (10) is normalized by N oc L D with 
L being the linear dimension of the system, A is propor- 
tional to (a/L) D , so that if a has similar values between 
2D and 3D we end up with a smaller A in 3D than that 
in 2D. So we can conclude that this is the main reason 
why the 2D case is more favorable than the 3D case. 

3.7 BCC lattice 

Let us turn to the BCC lattice. For BCC lattice, 
the density of states diverges around the center with 
D(E) ~ [ln(i?)] 2 , and the antiferromagnetic fluctuation 
is dominant near half-filling^ so we focus on the possi- 
bility of d-wave superconductivity. For the d-wave, we 
found that representation (~ xy etc, with the gap 
function A(k) ~ cos(k x + k y + k z ) — cos(k x + k y — k z ), 
etc) of Oh group has the largest Amux- 

In Fig. 17, we show Amox for the density varied over 
0.75 < n < 0.9 and f varied over -0.4 < t' < 0.4 for 
U — 4, 6, 8 with T = 0.03. Antiferromagnetic spin fluc- 
tuations are much stronger in BCC lattice than in the 
SC lattice as pointed out in ref.EF In fact, in this fig- 
ure, the truncated curves for AMax means that the sys- 
tem becomes antiferromagnetic (i.e., Uxo becomes unity) 
there. Hereafter, we focus on the case of U = 6, t' = 0.1, 
n = 0.75. 

In Fig. 18(a), we plot the momentum dependence of 
x(k, 0). We can see that the peak of %(k, 0) is located 
around fC-point (k = (tt, 7T, 7t)). In Fig. 18(b), we plot 
the frequency dependence of Imx(kMax, which shows 
that the frequency spread for the BCC lattice is similar 
to the case of the square lattice or the SC lattice. In 
Fig. 19, we plot AMax for c?-wave pairing as a function of 
T. AMax is again much smaller than that for the nearly 
half-filled square lattice, and is even smaller than that for 
the SC lattice. The reason for the former is mainly due 
to the (a/L) D factor discussed above, while the reason 
for the latter is because the maximum in t/ 2 |G| 2 in BCC 
lattice is smaller (see Fig. 20) than that in SC lattice 
(Fig. 16). 

Fig. 17. A plot similar to Fig. 6 for a BCC lattice, with the den- 
sity varied over 0.75 < n < 0.85 and t' varied over —0.4 < t' < 
0.4 for U = 4, 6, 8 with T = 0.03. Where the curves are trun- 
cated the system becomes antiferromagnetic (i.e., f/xirr becomes 
unity) . 



d- and p-wave superconductivity mediated by spin fluctuations 



Fig. 18. (a) 0) as a function of wave number for the Hubbard 
model on a BCC lattice with t' = 0.1, n = 0.75, T = 0.03 and 
U = 6. Note that if we take the energy dispersion (6), (it, it, it) 
and (n, 0,0) are equivalent. (b)Imx(kMax! ^) (normalized by its 
maximum value) as a function of uj/t for the Hubbard model on 
a BCC lattice for t' = 0.1, n = 0.75, (7 = 6 with T = 0.03 (left 
frame) or T = 0.04 (right). 



Fig. 19. The maximum eigenvalue of the Eliashberg equation 
(solid lines) and the reciprocal of the antiferromagnetic peak 
in x (dashed lines) against temperature for the Hubbard model 
on a BCC lattice with t' = -0.1, n = 0.75 and (7 = 6. 



Fig. 20. A plot for the squared Green function against k x and 
k y with k z = 0, 7r/2, 7r for the Hubbard model on a BCC lattice 
with t' = 0.1, n = 0.75, U = 6 and T = 0.03. 



3.8 FCC lattice 

We finally come to the FCC lattice. The density of 
states of FCC lattice diverges at the bottom of the band 
with D(E) ~ \hx{E + 4)| for t' = and - l/y/W+3 for 
t' = 0.50, so we can expect large ferromagnetic spin fluc- 
tuations for low densities of electrons. In fact, the pos- 
sibility of ferromagnetic gmtpd-state has been discussed 
using various approaches. O oeP According to our pre- 
vious study ,EF the ferromagnetic spin fluctuation is most 
dominant for n ~ 0.2 and t' ~ 0.5 in the weak coupling 
regime. Here, we focus on the possibility of p-wave pair- 
ing at low densities. In this case, we found that the case 
of the gap function A(k) ~ sin^ + k y ) has the largest 
Aiviax- (Since the Hubbard model has SU(2) symmetry, 
the T c does not depend on the direction of the d-vector, 
so that we may concentrate on the wave-number depen- 
dence.) 

In Fig. 21, we show the maximum eigenvalue of the 
Eliashberg equation for the density varied over 0.2 < n < 
0.4 and t' varied over 0.0 <t' < 0.6 for U = 2, 4, 6, 8 with 
T = 0.03. We can see that Aiviax is much smaller than 
unity even around t' ~ 0.5. 

To probe the origin of this behavior, we take the case 
of U = 2, t' = 0.5, and n — 0.3, which is also convenient 
becase the maximum value of x(k, 0) (see below) in this 
case takes a similar value as in the square lattice with tf — 
0.5 (Fig. 7(b)), which facilitates the comparison between 
the two cases. In Fig. 22(b), we show the momentum 
dependence of x(k, 0). We can see that the peak is indeed 
located around the ferromagnetic point (r). 

In Fig. 23, we plot AMax along with the reciprocal 
of the peak value of x(k, 0) as a function of T. AMax 
in this case is much smaller than that for the d-wave 
pairing in the SC lattice. This is again mainly due to 
the one-third V ^ . AMax is smaller even when compared 
with that for the p-wave pairing in the square lattice 
with t' = 0.5, although the maximum value of %(k, 0), 
the frequency spread of Imx(kMax. u>) (Fig. 22(a)), and 
£/ 2 |G| 2 (see Fig. 24) all take similar values between the 
two. In fact, the reason for this 2D-3D discrepancy can 
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again be traced back to the (a/L) D factor discussed in 
the previous subsection. 

Fig. 21. A plot similar to Fig. 6 for an FCC lattice, with the 
density varied over 0.2 < n < 0.4 and t' varied over 0.0 < t' < 0.6 
for U = 2, 4, 6, 8 with T = 0.03. Where the curves are truncated 
the system becomes ferromagnetic (i.e., J7xirr becomes unity). 



Fig. 22. (a)Imx(k[ V [ax) aJ ) (normalized by its maximum value) as 
a function of uj/t for the Hubbard model on an FCC lattice with 
t' = 0.5, n = 0.3, U = 2 and T = 0.03 (b)x(k, 0) as a function 
of wave number for the Hubbard model on an FCC lattice with 
t' = 0.5, n = 0.3, T = 0.03 and U = 2. Note that the T-point 
(0,0,0) and the X-point (tt,tt,tt) are equivalent if we take the 
dispersion (6). 



Fig. 23. The maximum eigenvalue of the Fliashberg equation 
(solid lines) and the reciprocal of the ferromagnetic peak of \ 
(dashed lines) in the Hubbard model on an FCC lattice with 
t' = 0.5, n = 0.3 and U = 2. 



Fig. 24. The squared Green function against k x and k y with k z = 
0, tt/2, tt for the Hubbard model on an FCC lattice with t' = 0.5, 
n = 0.3, U = 2 and T = 0.03. 



§4. Discussions and Summary 

To summarize, we have studied the possibility of spin- 
fluctuation mediated superconductivity in the single- 
band, repulsive Hubbard model for the d-wave channel 
on the (a) square lattice, (b) SC lattice, (c) BCC lat- 
tice, and for the p-wave channel on the (d) square lattice 
with large second-nearest neighbor hopping, (e) trian- 
gular lattice, and (vi) FCC lattice. We have shown that 
(i) d-wave instability mediated by antiferromagnetic spin 
fluctuations is stronger than p-wave instability mediated 
by ferromagnetic spin fluctuations both in 2D and 3D, 
and (ii)e?-wave instability in 2D is much stronger than 
that in 3D. 

We have given the physical reasons why the triplet p- 
pairing is unfavored as (i) the pairing interaction for 
triplet pairing is only 1/3 of those for the singlet pair- 
ing, (ii) the self energy correction is large (i.e., quansi- 
particles are short-lived). 

We have also traced back physical reasons why the 
superconducting instability in 3D is weaker than in 2D: 
if the momentum spread of x (that determines V^) in 
each momentum direction take similar values in 2D and 
3D, it makes the higher-dimensional 3D disadvantageous 
because of a structure of the Eliashberg equation. Thus, 
our conclusion is that, so far as the single-band Hub- 
bard model on ordinary lattices are concerned, d-wave 
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pairing in 2D square lattice is the "best" situation for 
the spin fluctuation mediated superconductivity, where 
T c can reach 0(Q.OQ1W) if we measure T c in units of the 
band width, W. In this sense, the layer-type cuprates do 
seem to hit upon the right situation. 

However, our conclusion has been obtained for the 
single-band Hubbard model. If wc look more extensively 
at 3D superconductors, the heavy fermion system, in 
which the pairing is also thought to be meditated by 
spin fluctuations, T c has a similar order of magnitude 
0(0.01 W ~ n.ooiWO, i-e., T C ~1K with W = a few 
hundred K.EP Since the present result indicates that 
T c in the 3D Hubbard model should be much smaller 
than that for the 2D Hubbard model, we can envisage 
that the heavy fermion system exploits a more favorable 
situation than in the single-band Hubbard model where 
the frequency and/or momentum spreads in x(k,w) are 
larger than those in the 3D Hubbard model. In fact, the 
standard models, such as the multiband periodic An- 
derson model, employed to describe the heavy fermion 
system are more complicated than the single-band Hub- 
bard model, and the frequency/momentum spreads in \ 
in such a model should be large enough to explain the 
superconductivity in the heavy fermion systems in the 
present context. It is an appealing future problem to 
explore how this can be so. 

§5. Acknowledgment 

We would like to thank Professor K. Ueda and Dr H. 
Kontani for illuminating discussions. R.A. would like 
to thank Dr S. Koikegami for discussions on the FLEX. 
R.A. is supported by JSPS, while K.K. acknowledges a 
Grant-in-Aid for Scientific Research from the Ministry 
of Education of Japan. Numerical calculations were per- 
formed at the Supercomputer Center, ISSP, University 
of Tokyo. 



[1] J. G. Bednorz and K. A. Miiller: Z. Phys. B 64 (1986) 189. 

[2] T. Moriya, Y. Takahashi and K. Ueda: J. Phys. Soc. Jpn, 59 
(1990) 2905; Physica C 185-189 (1991) 114. 

[3] T. Ueda, T. Moriya and Y. Takahashi: Electronic Proper- 
tics and Mechanisms of High-T c Superconductors, Tsukuba 
Symp., 1991, ed. T. Oguchi et al. (North Holland, Amster- 
dam, 1992) p. 145; J. Phys. Chem. Solids 53 (1992) 1515. 

[4] T. Moriya and K. Ueda: J. Phys. Soc. Jpn. 63 (1994) 1871. 

[5] P. Monthoux, A. V. Balatsky and D. Pines: Phys. Rev. B 46 
(1992) 14803. 

[6] P. Monthoux and D. Pines: Phys. Rev. B 47 (1993) 6069. 
[7] P. Monthoux and D. Pines: Phys. Rev. B 49 (1994) 4261. 
[8] R. Arita, S. Onoda, K. Kuroki and H. Aoki: to be published 

in J. Phys. Soc. Jpn. 
[9] N. E. Bickers, D. J. Scalapino and S. R. White: Phys. Rev. 

Lett. 62 (1989) 961. 
[10] N. E. Bickers and D. J. Scalapino: Ann. Phys. (N. Y.) 193 

(1989) 206. 

[11] T. Dahm and L. Tcwordt: Phys. Rev. B 52 (1995) 1297. 
[12] J.J. Deisz, D. W. Hess and J. W. Serene: Phys. Rev. Lett. 76 
(1996) 1312. 

[13] H. Kino and K. Kontani: J. Phys. Soc. Jpn. 67 (1998) 3691. 
[14] H. Kondo and T. Moriya: J. Phys. Soc. Jpn. 67 (1998) 3695. 
[15] J. Schmalian: Phys. Rev. Lett. 81 (1998) 4232. 
[16] M. Vojta and E. Daggoto: Phys. Rev. B. 59 (1999) R713. 
[17] T. Jujo, S. Koikegami and K. Yamada: J. Phys. Soc. Jpn. 68 
(1999) 1331. 



[18] T. Jujo and Yamada: J. Phys. Soc. Jpn. 68 (1999) 2198. 
[19] D. J. Scalapino, E. Loh, Jr. and J. E. Hirsh: Phys. Rev. B 34 
(1986) 8190. 

[20] S. Nakamura, T. Moriya and K. Ueda: J. Phys. Soc. Jpn 65 

(1996) 4026. 

[21] A. J. Leggett: Rev. Mod. Phys. 47 (1975) 331. 

[22] H. Tou, Y. Kitaoka, K. Asayama, K. Kimura, Y. Onuki, E. 

Yamamoto and K. Maezawa: Phys. Rev. Lett. 77 (1996) 1374; 

ibid, 80 (1998) 3129. 
[23] Y. Maeno, H. Hashimoto, K. Yoshida, S. Nishizaki, T. Fujita, 

J.G. Bednortz and F. Lichtenberg: Nature 372 (1994) 532. 
[24] T. M. Rice, M. Sigrist: J. Phys. Condens. Matter 7 (1995) 

L643. 

[25] W. Kohn and J.M. Luttinger: Phys. Rev. Lett. 20 (1965) 524. 
[26] M. Yu. Kagan and A. V. Chubukov: Pis'ma Zh. Eksp. Teor. 

Fiz. 47 (1988) 525. 
[27] A. V. Chubukov: Phys. Rev. B 48 (1993) 1097. 
[28] Y. Takada: Phys. Rev. B 47 (1993) 5202. 
[29] C. A. Kukkonen and A. W. Overhauser: Phys. Rev. B 20 

(1979) 550. 

[30] A. V. Chubukov and J. P. Lu: Phys. Rev. B 46 (1992) 11163. 

[31] R. Hlubina: Phys. Rev. B 59 (1999) 9600. 

[32] H. Takahashi [J. Phys. Soc. Jpn, 68 (1999) 194] on the other 
hand concludes that p-wave channel is most attractive for di- 
lute (n ~ 0.1) 2D Hubbard model on the t' = square lattice. 

[33] P. Monthoux and G. G. Lonzarich: Phys. Rev. B 59 (1999) 
14598. 

[34] R. Arita, K. Kuroki and H. Aoki: Phys. Rev. B 60 (1999) 
14585. 

[35] G. Baym and L.P. Kadanoff: Phys. Rev. 124 (1961) 287. 
[36] G. Baym: Phys. Rev. 127 (1962) 1391. 

[37] H. J. Vidberg and J. W. Serene: J. Low. Temp. Phys. 29 
(1977) 179. 

[38] N. D. Mermin and H. Wagner: Phys. Rev. Lett. 17 (1966) 
1133. 

[39] P.C. Hohenberg: Phys. Rev. 158 (1967) 383. 

[40] T. Koma and H. Tasaki: Phys. Rev. Lett. 68 (1992) 3248. 

[41] N. R. Werthamer: Superconductivity, ed. by R. D. Parks 
(Marcel Dekker, New York, 1969). 

[42] In liquid 3 He a p-wave superfluidity is mediated by ferro- 
magnetic spin fluctuations in an almost ferromagnetic system, 
which is usually interpreted in terms of the paramagnon ex- 
change for a system of particles with hard-core interactions. 
Two-dimensional liquid 3 He has also been studied, but the 
paramagnon mode is shown to be very sloppy [M. Ogura and 
H. Namaizawa, J. Phys. Soc. Jpn 66 (1997) 3706]. While the 
fermions move in a continuous space in liquid 3 He, here the 
ferromagnetic fluctuation becomes dominant because the den- 
sity of states can be large depending on the lattice structure, 
which is an effect specific to lattice fermions. 

[43] R. Hlubina, S. Sorella and F. Guinea: Phys. Rev. Lett. 78 

(1997) 1343. 

[44] T. Hanisch, G.S. Uhrig, E. Miiller-Hartmann: Phys. Rev. B 

56 (1997) 13960. 
[45] M. Sigrist and K. Ueda: Rev. Mod. Phys. 63 (1991) 239. 
[46] M. Ulmke: Euro. Phys. J. Bl (1998) 301. 



This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig2.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig3.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig4.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig5.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig6.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig7.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig8.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig9.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "figl0.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "figll.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "figl2.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "figl3.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "figl4.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "figl5.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "figl6.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "figl7.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "figl8.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "figl9.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig20.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig21.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig22.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig23.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig24.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig25.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig26.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig27.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig28.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig29.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



This figure "fig30.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0002440v 1 



